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L.  INTRODUCTION 

Evaluating  the  role  of  electric  currents  in  the  energy  storage  process 
within  coronal  and  chromospheric  magnetic  fields  has  been  limited  mainly 
to  comparing  field  lines  extrapolated  from  observed  longitudinal  photo- 
spheric  values  to  emission  features  assumed  to  be  field-aligned  (  Tanaka  and 
Nakagawa,  1973;  Hagyard,  1976;  Levine,  1976;  Krall,  et ,  al.,  1978). 

These  calculations  have  been  done  using  arbitrary  and  overly  simplistic 
constant-a  (  =  p^J/B,  where  J  is  the  current  density,  _B  the  magnetic 
field,  and  Uq  the  permeability  in  free  space)  assumptions  concerning 
current  distributions.  Measuring  the  transverse  component  of  the  photo- 
spheric  field  allows  a  unique  opportunity  to  obtain  information  about 
a-distribution  and  thus  the  Lorentz  forces  exerted  on  the  overlying  layers. 
Unfortunately,  utilizing  an  observationally  deduced  a-distribution  to 
compute  the  force-free  field  in  the  region  above  the  observing  plane  is 
in  general  an  intractable,  nonlinear  boundary  value  problem.  However, 
a  method  exists  for  determining  net  Lorentz  forces  and  magnetic  energy 
within  a  force-free  region  without  explicit  solution  of  the  field  equations 
(eg.  Molodensky,  1974).  The  purpose  of  this  paper  is  to 
examine  the  utility  of  such  a  method,  within  the  limitations  of  available 
vector  magnetic  field  observations,  in  determining  energy  exchange  between 
photospheric  fields  and  the  overlying  atmosphere.  Such  a  method  may  be 
potentially  useful  in  real  time  for  identifying  modes  and  levels  of  energy 
buildup  in  the  preflare  active  region. 

A  discussion  of  pertinent  equations  is  presented  in  Section  2.  In 
order  to  properly  evaluate  the  results,  measurement  and  calibration  un¬ 
certainties  are  discussed  in  sections  3  and  4.  Although  the  interpretation 


of  transverse  measurements  is  subject  to  significant  uncertainty, 
correlative  observations  may  be  utilized  to  test  the  integrity  and 
resultant  conclusions  of  the  measurements  (e.g.,  photospheric  proper 
motion  may  be  found  to  be  consistent  with  generation  of  the  J  x  B  forces 
inferred  from  the  measured  fields).  Section  5  presents  application  of  the 
previously  discussed  formalism  to  two  sets  of  Marshall  Space  Flight  Center 
Vector  Magnetograph  data,  a  flare  productive  and  migrating  delta-configura¬ 
tion  spot  on  6  April  1980  and  a  large  but  rather  simple  (beta-gamma)  region 
on  15-16  September  1980. 
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2.  Theoretical  Formulation 

The  vector  sum  of  all  electromagnetic  forces  acting  within  a  volume 
V  above  the  photosphere  can  be  written  as: 


F 


1 

c 


J  x  B  dV 
■*V 


(2-1) 


where  c  is  the  velocity  of  light  and  we  assume  electric  fields  to  play 
a  negligible  role  due  to  the  high  electrical  conductivity.  Using  Max¬ 
well's  equations 


J 


c 

4tt 


V  x  B 


(2-2) 


and  V  •  B  =  0  , 


(2-3) 


Equation  (2-1)  can  be  written 


(V  x 


B)  x  B  +  B  (V  •  B)  dV 


) 


(B  •  V)B  -  |VB2  +  B(V  •  B)  dV 


(2-4) 


The  integral  in  Equation  (2-4)  can  be  identified  as  the  divergence  of 
the  second  rank  Maxwell's  stress  tension 


T  ■  ^  (I  •  B 


B2) 


(2-5) 


whose  elements  are 
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7^  (b  B  B2 

4ir  (  i  j  2  ij 


(2-6) 


where 


6ij  -1  i‘j 

-  0  i  +  j  . 


(2-7) 
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Thus,  the  volume  Integral  Equation  (2-4)  can  be  replaced  by  an  integral 
over  the  surface  S  containing  V: 


F  -  I  7  *  T  dV 
JV 


•  T  dS 


(2-8) 


where  n  is  the  outward  direction  unit  normal  to  S.  In  a  similar  fashion, 
the  net  electromagnetic  torque  X  acting  within  a  volume  can  be  shown  to 


~  r  x  (J  x  B)  dV 

Jv 


i  x  (V  •  ij  dV 
■*V 

n  •  (T  x  r)  dS 
■  s 


where  “  ^ji  use<*  to  arrive  at  Equation  (2-9)  and  r  is  the  position 
vector  of  a  lever  arm.  Thus,  if  a  volume  enclosing  chromospheric  coronal 
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fields  is  appropriately  chosen  to  give  zero  contribution  to  Equations  (2-8) 
or  (2-9)  on  all  but  Lire  photospheric  surface,  a  complete  specification  of  (1 
on  that  level  is  sufficient  to  determine  a  lower  limit  (net)  to  the  magnetic 
forces  or  torques  exerted  on  the  volume.  The  work  done  by  the 
photospheric  fields  and  any  resulting  magnetic  energy  buildup  cannot,  of 
course,  be  determined  by  this  method.  In  the  absense  of  chromospheric 
or  coronal  magnetic  field  data,  these  can  only  be  inferred  indirectly 
(e.g.  by  changing  fibril  structure  or  energy  release). 

The  net  forces  and  torques,  and  total  magnetic  energy  (which  is 
derived  similarly)  can  be  expressed  in  Cartesian  coordinates  using 
longitudinal  (z  direction)  and  transverse  (x-y  plane)  field  measurements. 
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3.  Measurements 


The  MSFC  Vector  Magnetograph  utilizes  a  .1258  Zeiss  filter,  tunable  to 
10  m&  increments  within  the  FeT  5250.28  line.  All  observations  discussed  here 
are  centered  at  X  *  60,  90  or  120  m8  in  the  blue  wing  of  that  line.  The 
measured  components  of  the  magnetic  field  are  assumed  to  be 


B  =  k  C. (AX)  P  (AX) 

z  1  y 


2  2 

B  +  B 
x  y) 


k2  C2(AA)  Pq (AX) 


(3-1) 


(3-2) 


where  the  fractional  circularly-polarized  intensity 

Pv  “  K  '  h}/  (xi  +  h]  • 


and  linearly  polarized  intensity 


'  2 

U  +  R  1 


and 


(3-3) 


(3-4) 


X3-  *4 


)/( 


X3  +  I4 


R  = 


X5  ~  *6 


I.  +  1,1 

5  6  j 


(3-5) 


are  functions  of  the  MSFC  Vector  Magnetograph  raw  intensity  counts 
described  in  Table  I.  The  direction  4)  of  the  tranverse  field  where 


i 


and 


B  ■  B„  cos  d> 
xT 


B  =  B  sin  4>, 


0-6  ) 


is  found  as 


4>  =  .5  tan  ^  (U/R).  (3-7) 

The  values  for  were  determined  from  theoretical  calibration  curves 

derived  from  the  kjeldseth  Moe  solutions  to  the  radiative  transfer  equations 
(Kjeldseth  Moe,  1968)  for  penumbral  and  umbral  atmospheric  models. 

The  parameter  k  represents  a  scaling  factor  which  adjusts  measured 
polarization  for  loss  of  contrast  and  was  determined  by  scaling  the  maximum 
B  magnitude  with  that  of  Mt  Wilson  measurements.  Corrections  for  linear 
and  circular  crosstalk  in  the  MSFC  data  are  discussed  in  the  following 
paragraphs. 

Average  background  (quiet  sun)  measurements  of  Py  and  are  normally 
used  to  look  for  linear  and  circular  polarization  biases.  Typically,  no 
circular  bias  is  present;  however,  random  linear  polarization  biases  are 
found  and  are  attributed  to  the  magnetograph's  correlation  tracker  subsystem. 
Figure  1  a-d  shows  an  example  of  the  distribution  of  weak  field  polarization 
values  used  in  determining  these  biases.  Figure  2  a-d  shows  corresponding 
distributions  after  correcting  for  bias. 

Circular  crosstalk,  the  effect  of  coupling  circular  into  linear 
polarization  measurements  in  strong  field  areas,  must  also  be  considered 
in  the  calibration  process  to  eliminate  biasing  <J>.  Figure  3  a,b  are 
scatter  plots  of  U  and  R  vs  P^ ,  with  the  linear  least  squares  fits  showing 
the  artificial  correlation  induced  by  crosstalk.  A  correction  factor  has 
been  introduced  in  Figure  4  a,b  to  eliminate  the  correlation. 
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4.  Error  Analysis 

The  effects  of  signal  to  noise  and  calibration  must  be  estimated  in 
order  to  determine  the  integrity  of  the  measurements.  This  section 
describes  the  treatment  of  noise  in  the  linear  arid  circular  polarization 
measurements  used  to  derive  longitudinal  and  transverse  fields. 

The  physical  parameters  to  be  evaluated  are  Equations  (2-10)  through 
(2-13)  and  direction  of  horizontal  force: 


0  *=  tan 


-1 


F 

JL 

F 

x 


Gaussian  statistics  are  used  to  determine  the  standard  deviation  at  each 
pixel  i, j : 


Ay  B  -  Ax  B  )  }  AS2,  (4-1) 
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B  o^  +  B  0*  +  B  at 
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1 

2 

1 

F 

rr  O  i  jf  /r  O 

1  +  (F  /F  )2 

L  y  x  J 

F2 

L  x 

+  , —  o_z 

F  r4  F 

y  F  x 

1  X 

(4-3) 


(4-4) 


Equations  3-1  through  3-9  are  used  to  determine  the  field  uncertainties 


on  MSFC  data  to  calculate  Equations  2-10  through  2-13  and  8,  as  well  as 
standard  deviations  4-1  through  4-4.  Input  from  data  files  are  raw 
intensitites  1^  through  1^  and  processed  fields  Bz»  B^,  B^  (i.e. ,  after 
calibrations  and  elimination  of  the  180°  ambiguity  in  transverse  field 
direction) . 

The  standard  deviation  in  the  raw  data  (aT  )  produces  the  spread  in 

Ai 

the  zero-field  distributions  (a  and  0  )  in  Figure  2.  Equations  4-5 

PQ  PV 

through  4-7  can  therefore  also  be  written  in  terms  of  measured  OP^  and 
Op  ,  and  are  written  so  within  the  present  code. 

*  *7 


5.  Results 


In  an  earlier  calculation  (Krall,  et  al.  1982),  a  small,  migrating 
and  flare  producing  portion  of  Boulder  SESC  AR2372  (Figure  5)  was  eval¬ 
uated  using  the  equations  of  Section  2.' 

For  the  48  hour  period  containing  the  maximum  photospheric  proper  motion 
and  flare  activity.  Figure  6  shows  the  direction  of  horizontal  magnetic 
force  for  the  area  containing  the  small  westward  moving  positive  spot. 

The  calculated  force  was  directed  within  30  degress  of  the  observed  motion 
of  the  spot  during  the  time  of  maximum  spot  motion  and  flare  activity, 
changing  significantly  (and  decreasing  somewhat)  outside  of  that  time  period. 
This  is  as  expected  if  the  photospheric  motions  act  to  magnetically  transfer 
forces  to  the  overlying  atmosphere.  Such  forces  can  of  course  be  dissipated 
rapidly  or  stored  in  a  higher  energy  magnetic  configuration.  The  later 
possibility  is  more  likely  here  due  to  the  high  flare  activity  of  the  period. 

Although  AR2372  was  an  extremely  interesting  region  from  the  viewpoint 
of  tii is  analysis,  several  factors  combined  to  complicate  the  interpretation: 

1.  There  are  possibly  strong  field  saturation  effects. 

2.  Region  complexity  made  difficult  the  task  of  elim¬ 
inating  the  180  degree  ambiguity  in  the  direction 
of  the  tranverse  field  component. 

3.  Inadequate  observations  at  greater  than  60  mX  from 
the  line  center  increased  the  possibility  of 
magneto-optical  effects  (West  and  Hagyard,  1983). 

4.  There  was  insufficient  time  resolution  for  r>  i 
analysis  of  magnetic  energy  buildup  or  dissipa¬ 
tion  for  such  a  rapidly  evolving  region  (During 
the  period  an  M  or  X  class  flare  occurred  on  the 

,  average  every  5  hours  and  C  class  flare  every 

I  8  hours). 
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It  was  therefore  decided  that  a  region  of  less  spatial  and  temporal 
complexity  was  needed  to  properly  evaluate  the  potential  of  this  type 
of  analysis.  Boulder  SGSC  If  AR2665  (Figure  7),  a  simple  large  spot  region 
on  its  second  disk  transit,  was  well  observed  at  90  and  120  m&  on  Sept¬ 
ember  IS  and  16,  had  little  indication  of  saturation  at  those  wavelengths, 
and  was  by  its  single  spot  nature,  relatively  easy  to  treat  for  its 
transverse  field  direction  ambiguity.  Its  most  significant  characteristics 
were  a  large  leader  spot  (of  positive  polarity)  and  a  lack  of  pronounced 
magnetic  complexity. 

During  the  first  few  days  of  its  rotation  across  the  visible  disk, 

AR2665  did  exhibit  some  complexity  and  some  evidence  of  sheared  fields, 
and  produced  several  relatively  small  flares  near  the  weakly  sheared 
neutral  line.  Flare  history  (Figure  8)  and  magnetic  history  (Figure  9) 
both  reached  maxima  from  September  11  -  13. 

The  region  had  been  simplifying,  and  continued  to  simplify  during  the 
period  of  interest.  Between  15/1349UT  and  16/1359111'  (Figure  7a  and  c)  the  BL 
gradients  along  the  neutral  line  had  decreased,  the  field  strengths 
(both  and  B^)  had  declined  near  the  neutral  line,  and  the  B^,  config¬ 
uration  had  become  more  potential  appearing  (orthogonal  to  the  B^  =  0  line) 
in  orientation. 

The  more  frequent  flare  activity  occurred  before  the  times  of  the 
Lorentz  force  evaluation  (Figure  8),  with  only  one  Cl  and  2  subflares 
without  associated  x-rays  occurring  during  the  72  hours  after  the  beginning 
of  the  analysis.  AR2665  had  been  reported  as  a  beta-gamma  region  (somewhat 
magnetically  complex)  but  by  the  15th,  observations  of  the  magnetic  field 
indicated  it  to  be  a  simple  beta  (simple  bipole)  magnetic  configuration. 


Figure  10  shows  physical  parameters  (Equations  (2-11]  through  (2-13) 
tor  a  2’ i  hour  period  during  the  decay  of  region  AR  2665.  Points  to  be 


emphasized  are: 


1.  Signal  to  noise  is  seen  in  all  cases  to  be  small 
compared  to  time  variation  and  observing  wave¬ 
length  variation. 

2.  Although  wavelength  variation  is  only  available 
at  one  time  point  in  Figure  10  (actually  a  twenty 
minute  period  at  2000  UT  on  15  September),  the 
results  are  consistent  with  a  measurement  de¬ 
gradation  as  one  approaches  line  center,  e.g., 

a  negative  vertical  force  at  60  raft  requires 
unlikely  submerging  magnetic  flux  or  upward 
counter-forces.  Such  degradation  is  consistent 
with  magneto-optical  effects;  in  fact,  detailed 
analysis  for  this  region  shows  that  reliable 
transverse  field  directions  may  be  obtainable 
only  at  AA  120  m&  (West  and  Hagyard,  1983). 

In  addition,  the  tightening  spiral  configuration 
suggested  by  the  divergence  in  torque  with 
decreasing  observation  wavelength  is  also 
consistent  with  magneto  optical  effects. 

3.  In  the  declining  phase  of  AR  2665,  available 
energy  is  seen  to  decay  by  10^2  erg s/day.  If 
this  is  not  an  observational  effect,  significant 
energy  release  over  a  relatively  short  time 
period  (10^~  ergs/ few  hours)  might  be  claimed 

t:o  have  been  observed. 

4.  Upward  vertical  force  decreases  during  the  24 
hours;  this  could  be  interpreted  as  a  declining 
rate  of  emerging  flux.  A  declining,  yet  always 
positive  rate  of  emerging  flux  is  to  be  expected 
during  the  demise  of  an  active  region. 

(This  upward  energy  flux  would  be  more  than 
compensated  for  by  dissipative  processes  during 
the  decay  phase,  whereas  the  converse  would 
be  true  during  region  growth.) 


6.  Concluding  Remarks 


We  have  developed  and  utilized  a  procedure  for  evaluating  the 
potential  role  of  3-component  photospheric  magnetic  field  measurements 
for  studying  the  magnetic  energy  buildup  process  within  solar  active 
regions.  We  have  applied  the  equations  to  two  regions  and  the  results 
appear  encouragingly  consistent  with  expectations,  both  as  a  function 
of  time  (physical  variations)  and  wavelength  (observational  variations). 
Further  evidence  is  presented  that  these  type  observations  must  be  made 
at  or  greater  than  90  n&  from  line  center. 

Important  questions  as  yet  to  be  answered  during  the  validation  pro 
cess  are: 

1.  Does  the  measured  vertical  magnetic  force 
increase  during  early  active  region  growth 
(when  emerging  flux  is  most  apparent), 
reach  a  maximum,  and  then  decay  within 
the  proper  time  frame?  This  is  probably 
not  answerable  within  the  lifetime  of 
single  active  regions,  but  can  be  approached 
by  studying  regions  in  each  phase  of  growth. 

2.  Do  proper  spot  motions  (or  rotations)  show 
up  as  increased  horizontal  magnetic  force 
(torque),  and  can  these  motions  now  be  more 
quantitatively  linked  to  magnetic  energy 
buildup  and  release? 


I  7 


C  THIS  PROGRAM  CALCULATES  the  3  COMPONENTS  OF  NET  FORCE,  TOROUE, 

C  ENERGY  AND  THE  SIGNAL/NOISE  UNCERTAINTIES  ASSOCIATED  UITH  THEM 

C  INPUT  FROM  TAPE  9  IS  THE  OBSERVED  PHOTOSPHERIC  BX,PY  AND  BZ 

C  INPUT  FROM  TAPE  10  IS  11,1?,. ..16  RAW  DATA  CORRESPONDING  TO  FIELD 

C  STRENGTHS  FROM  TAPE  9 

DIMENSION  S16T(41,41),SIGE<41  ,41  ),P(3,*1  ,41)  ,IR  (6)  ,AM  <10> 
DIMENSION  F (3) ,YORIGN(10) »XORIGN(10) 

DIMENSION  C10LD (6) ,C1NEW(6) ,C20LD(6) ,C?NEW<£ ) 

DIMENSION  P  TPH I (41  ,41) 

REAL  KN  EW ( £  ) 

REAL  KOL  D ( 6 ) 

REAL  IT(6 ,41  ,41) 

REAL  KFACT 

REAL  I56,I34,I56P,I34P,I12,I12P 
C 
C 

DATA  YORIGN/19. ,16., 16. , 1 9. , 1 5 . , 5* .0/ 

DATA  XORIGN/25.,23.,24.,23. ,25.,5*C./ 

C 

C 

DATA  IR/1 ,2, 3,4, 5,0/ 

C  UHBRAL  MODEL 

DATA  Cl  NEW/4  374. ,2707. ,31 95., 4374. ,4  374.  ,Q./ 

DATA  C 1 0LD /6690., 3 182., 43 11., 6690., 6690. ,0./ 

DATA  C?NEW/362£., 2568., 2910. ,3628., 362$. ,0./ 

DATA  C20LD/4437., 2591., 3191., 4437. ,4437.,u./ 

DATA  KNEW  /7. 18, 9. 16, 8. 56, 7. 03,6.72,0./ 

DATA  KOLD  /4 . 7, 7.8 ,6 . 35 , 4 .6 , 4 . 4 , C. / 

DATA  AM/1 • ,Q»,C.,1 . , 0 • , C  .  ,4*0./ 

DS=2.5*725E5 
PI=3. 14159 
DTOR*PI/ 180. 

C  MAX  PI/2  UNCERTAINTY  IN  DIRECTION 

SPHIMX=PI/2. 

FRPI=4.*PI 

C  DEFINE  SIGPO,  SET  INITIAL  GUESS  ON  SIG3**2 

SIGPQ=.001 
SIG32=17C. 

C 

C 

DO  1C  KTAPE1=1,6. 

C  1 

c 

C  INITIAL  TIKE  FOR  FINDING  INTENSITY  VA R I ANC E S  , S UC E E D I NG  TIKES  FOR 

C  CALC  PHYSICAL  SIGMAS 

KTAPE=KT4PE1 

IFUTAPE1  .GT.1)  KTAPE  =  KTAPE1-1 
IF<KTAPE.EQ.1>  REWIND9 
C  RESET  VARIABLES 

GO  TO  10? 

109  CONTINUE 
C  READ  DATA 

GO  TO  120 
129  CONTINUE 

C  LOOP  THROUGH  DATA  ARRAY,  CALCULATING  UNCERTAINTIES 

GO  TO  140 
149  CONTINUE 
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C 

DT0FY=F(2)  /  <F(?>*«2  ♦  F(3>**2) 

?TDFX=-F (3)  /  (F(2>**2  ♦  F(3>**2> 

TAUTOT=TAU/TAUTOT 
E  TOT=£  N /E  TOT 
FXTOT=F (2) /FXTOT 
FYTOT=F<?) /FYTOT 
TAU=TAU*DS*D$*DS/FRPI 
EN=EN*DS*OS*DS/FRPI 
DO  2  1=1,3 

2  F (1 )=F (I)*DS*DS/FRPI 

FHMAG=SQRT(F(2)*F<2>  ♦  F(3)*F(3>) 

C  ROTATE  BACK  TO  X-Y  COOROS 

C  FX=F(2>*COS(2O.*DTOR)-F(3)*SIN<20.*DTOR) 

C  F(3)  =  F(2)*SIN<20.*DTOR)-fF<3)*COS<20.*DTOR> 

C  F(2)=FX 

THETA=F(3)/F(2) 

THET A=A TAN (THETA I/DTOR 
IF(ITP.EQ.O)  GO  TO  9 

SIGTH2=SIGFX*DTDFX**2  ♦  S IG F Y*DTD F Y**2 
SIGDEG=SQRT(SIGTH2)/DTOR 

SIG32P*  S1GPO  *  S1GPQ  *  (ISTOP-ISTART+1 >  *  ( JSTOP-JSTART+1 ) /DPQDI2 

IF (KTAPEl .EQ.1)  SIG32=SIG32P 

SlGTdO,  JO=SQRT(SIGT(IO,JO))*OS*DS 

S1GE(1C, JO)*SQRT<SIGE<IO»JO))*OS*OS*OS/FRPI 

SI6FZ*SQRT(SIGFZ)*DS*DS/FRPI 

S1GFX=SQRT(SIGFX)*DS*DS/FRPI 

SIGFY=SQRT(SIGFY)*OS*OS/FRPI 

9  CONTINUE 
WRITE(6,6) 

6  FORMAT (/// ,SOH  TORQUE  ENERGY  FZ  FX  FV  ) 

WRITF (6,7)  TAU,EN, <F (I) , 1*1  ,3) 

WRITE (6 ,7)  SIGT (I G  »  J  0 ) ,  S  I  GE (1 0 , J  0 ) *  SIGFX,SIGFY t  S I G  F  2 
WRITE(6,7)  TAUTOT.ETOT, FXTOT, FYTOT  , FHM  AG 
WRITE (6 ,8)  SIGDEG,S1G32,SIG32P,THETA,SIGPQ 
8  FORMAT (// 

%  1 4H  sigma  THETA*  ,E10.3  »25H  D  SIGMA  INTENSITY  USED,E10.3, 

S25H  SIGM  INTENSITY  CALC*, El C.3 , 1  OH  THETA* , FS •  1  , 

S12H  SIGM  PQ=  , E8  .2  ) 

7  FORMATE5E1C.3) 

10  CONTINUE 
STOP 

C  *  *  *  *  BLOCK  IOC  *  *  *  * 

C  (BLKC) 

C  INITIALIZE  VARIABLES 

C 

100  CONTINUE 

C  RESET  VARIABLES 

ITP=IR(KTAPE) 

WRITE(6,1Ct) 

103  FORMAT (1M1  ) 

WRITE(6,1C2>  KT  APE ,  ITP 
1 02  F0RMAT(2I1C) 

SIG1?*SIG32 
S IGZ2*S  IG32 
SIG42=SIG32 
SIG5?=SIG3? 


7  vr*  y 
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SIG62=SIG32 

C  I KNOWN*0  FOR  K  ONLY  KNOWN',  1  FOR  K,C1,C2  KNOWN 

IKNOWN=1 

KFACT=KNEU(KTAPE> 

C  K  F  A  CT  =  4  . 

C  C1=100C. 

C  Ct=1 COO. 

C1*ClNEW<KTAPE> 

C2*C2NEW<KTAP£> 

JG=X0RIGN<KTAPE> 

lC*YORIGN(KTAPE) 

TAU=0. 

£N=0. 

S1GTCIC, J0)=0. 

SIGEdC,  JO>=0. 

SIGFX*C’. 

SIGFY=G. 

SIGFZ=C. 

SIG02*C. 

DPGDI2-0. 

DO  101  1*1,3 

101  F ( I ) *0 • 

THETAsO. 

S1GTH2=C. 

ETOT-O • 

T  AUT  OT -0 • 

FXT0T=C. 

FYT0T=0.0 
GO  TO  109 

C  *  *  *  *  BLOCK  120  *  *  *  * 

C  (BLKC) 

120  CONTINUE 
C 

C  READ  DATA  FILES 

READ (9)  P , N , C ( (B (K ,1 ,  J  ) ,  J  *1 , N ) ,1*1 »M  ) *K  »1 ,3) 

C 
C 
C 

WRITE(6,124) 

124  FORMAT (1  Hi ) 

WRITE <6 ,123>(CB<1, I,J),J*1,39,2), 1*1 ,39,2) 

123  FORMAT (2CF6.0) 

DO  125  1*1  ,  M 
DO  125  J * 1  , N 

eTPHI Cl ,J)=ATAN<BC3,I,J)/p(2,I,J))/DT0R 

IF(B(3,I»J>.LT.0..AND.B(2,I  ,J).LT.O.)  BTPH I ( I , J ) *BTPH I ( I , J ) -IE 
IF(B(3,I,J).GT.0..AND.9(2,I  ,J).LT.O.)  BTPHI(I,J)*BTPHI(I,J)*1S 

125  CONTINUE 
WRITE(6,12A) 

WRITE(6,123)(<BTPHI(I,J),J=1,39,2),I=1,39,2) 

WRITE  (6,124) 

C 

c 
c 

PEWINO  10 

IFUTP.fcQ.C)  GO  TO  122 
DO  121  I T  A°E  *1 , 1  TP 


C->  Cl 
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IF  (AM(ITP).EQ.I.)  READ  C1D((  (IT  (K,I#J),I=H,  1,-D  t  J=N,1,-1)  »K-1,6) 
IF(AM(ITP).EQ.O.)  *EADCl0)CCCIT<K,I#J)tI»1 ,M),J=1 ,N),K=1,6> 

121  CONTINUE 
122  CONTINUE 
GO  TO  12? 

r  *  *  *  *  BLOCK  14  0  *  *  *  * 

J  (BLK  0 ) 

140  CONTINUE 

C  LOOP  THROUGH  DATA  ARRAY 

C 

ISTOP=P 
JSTART=1 
JSTOP=N 
ISTAPT=1 
ICHK=I START+3 
JCHK=J  START*? 

C  RECOVER  RAW  DATA 

DO  141  1=1 , M 
DC  141  J  =1 , N 
DO  141  K  C  AL  =  1 ,6 

141  IF  (IT(<CAL,I,J).LT.O.)  I T ( K CAL , I , J > =IT C KC AL , I , J > *65 536 . 


USE  CORRECT  OR  UPDATED  CONSTANTS 


C0R1=1 . 
C0P2=1 . 


C0R2=SQRT(tCNEW  (I TP)  /KOLD  (ITP)  )*C2NEW  ClTP)  /C20LD  (I TP) 

COR1=KNEW(ITP)*ClNEW<lTP)/<KOLD(ITP)*ClOlD<ITP>) 


DO  148  I=ISTART,ISTOP 
D  I  =  I- 1 C 

DO  147  J=JSTART, JSTOP 
D J=J-J  3 

0(1  , 1, J)=B<1, I, J)*C0R1 
P<2,I,J)*B<2,I,J>*COR2 
B<7, 1, J)=B(3,I,J>*C0R2 
3X  =  C<2,I,J) 

B  Y  =  B (2 , 1  ,  J  ) 

BZ=B  <1 , I  ,  J  ) 


CALCULATE  PHYSICALS 
GO  TO  1  GOC 
1009  CONTINUE 

CORRECT  SIGNS  OF  8  IF  N£CES$ARY(  Y  UP,BX  TO  RIGHTfBZ  OUT  OF  PLANE) 
BX  =  -(P  <2, I , J>*COS(20.*DT0R)*B  <3,1 ,J  >*SIN  C20.*DT0R) > 
r  Y  =  -  8  <2 , 1 ,j)*SIN(20,*DTOR)*B (3 , I , J )*COS < 20.  *DTOP> 

B  I  *B  (  1  »  I  ,  J  ) 

PT=SGRT (GX*8X*RY*BY) 

IF(ITr.EO.C)  GO  TO  147 


1 
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C  IF  CALIBRATION  CONSTANTS  NOT  KNOWN,  CALCULATE  THE"'! 

C 

IF ( IKNOWN . EQ.O.AND . I .EQ .1  START . AND . J .EC  .  J ST  APT  )  GO  TO  101D 
1 01 0  CONTINUE 

C  CALCULATE  INTENSITY  DIFFERENCES  ,SU!*S 

156=  IT(5,I,J)-IT(6,I,J) 

I56P»IT(5,I,J)  +  IT(6,  I,J) 

134=  1T<3, I,J)-IT(4,I, J) 

I34P  =  IT(3,I,J)+IT(4,1,J) 

112=  IT(1,I, J)-IT(2,I*J) 

I12P*  ITC1 ,1, J)4ITC2,I, J) 

C  CALCULATE  U , V , S IGU , S IGV 

U=  1?4 / 1 34  p 
V=  I56/I56P 
UVMAG=U*U+V*V 

U2V2*  (U/l34p>**2  ♦  <V/I56P)«*2 
C  DPCDI2=DPQDI2  ♦  2. *KFACT*KFACT*U2V2  /  (U*U+V*V> 

DPQDI2=DPQDl2  ♦  2.  *  U2V2  /  (U*U+V*V) 

$IGU=  <1.-U)**2  *  CSI632  ♦  SIG42)/ (I34P*I34P) 

SIGV*  (1.-V)**2  *  (SIG52  ♦  SIG62)/(I56P*I56P) 

C  IF(KTAPE.EQ.4) 

C  S WRITE (6, 144)  156 ,1 56P, 134 , 1 34P , 1 1 2 , 1 1 2P , U , V ,UVK AG ,U2V2  ,S IGU  ,  S I GV 
144  FORMAT  (12E10.3) 

C 

C  DETERMINE  DIRECTION  PARAMETERS 

GO  TO  102C 
1029  CONTINUE 
C  GO  TO  1050 
C1C59  CONTINUE 
C  CALCULATE  B  SIGMAS 

GO  TO  1030 
1039  CONTINUE 

C  CALCULATE  PARAMETER  SIGMAS,  CHECK  FOR  LARGE  TERMS 

GO  TO  1040 
1049  CONTINUE 
147  CONTINUE 
14?  CONTINUE 
GO  TO  149 

C  *  *  *  *  BLOCK  1000  *  *  *  * 

C  (BLK  140 

1C00  CONTINUE 
C 

C  CALCULATE  FORCES, ETC. 

C  TORQUE  =  ez  *  (BX  *  DY  -PY*DX) 

T=BZ*(BX*DI-BY*DJ) 

E=BZ*  (BX * D  J  ♦  B  Y*D I  ) 

F1=BZ*BZ-BX*BX-BY*9Y 

F2=BZ*BX 

F3=PZ*PY 

C  SUM  OVER  ARRAY 

EN=EN+E 
TAU=TAU*T 
F<1)=F(1)+F1 
F  (2  )  =  F  (2  )  ♦  F2 
F(3)*F<3>*F3 
ETOT  =  ETOT-»  ABS  (E) 

TAUT0T=TAUT0T-»A8F  <T) 


FXT0T=FXT0T4ABS ( F 2 ) 

FVTOT  =  F YT0T4ABS  <  F  3  > 
rc  TO  1 009 

C  *  *  +  *  BLOCK  1010  *  *  *  * 

C  (BLK  140) 

101C  CONTINUE 
I  FLAG  =  C 
V  AR 1 =0 . 

VAR2=0 • 

1015  CONTINUE 
C 

C  CALCULATE  Cl, C2, GIVEN  KFACT 

I C  OUNT  =  0 
JCOUNT=G 
RATI01=C. 

R  AT 1 02  =0 • 

BLIK-500. 

PLIM  =  10  . 

DO  101S  110=1 ,K 
DO  1017  J 1 0  =  1 f  N 

Fzc=e<i ,iio,jio) 

IF(eZC.LT.BLIM.AND.BZC.GT.-BLIM)  GO  TO  1C11 
BZP  =  KFACT  * 

1  C I T  C 1 . 110, J10)-IT<2, 110, J1C) )/<ITC1 ,I10,J1C)MTC2,I10, J10>> 
RAT1=B7C/BZP 

IF(IFLAG.EQ.I)  V AR 1 =V AR 1 ♦( R AT 1 - C 1 ) * *2 
C  WRITE  (6  ,9  011  )I10,J10,RATl  ,BZC  ,BZP,IT<1,I10,Jl0),IT(2,H0,J1G) 

C9011  FORMATC2I5.7F10.3) 

RATI01=RATI0l*BZC/BZP 
ICOUNT=I COUNT+1 
1  C 1 1  CONTINUE 

BTC=SQRT(B(2,I10,Jl0>**2  4B<3tI10,Jl0>W> 

I  F  (BTC  .LT.BLIM.AND .BTC .GT.-BLIM)  GO  TO  1012 

U=(IT(?,I10,J10)-IT(4,I10,J10))/(IT(3,I10,J10)*ITC4,I10,J10)) 

V=(IT(?,I1C,J1C)-IT(6,I10»J1Q))/(IT(5,I10,J10)4IT(6»I1C,J10)) 

UVMAG=U*U*V*V 

BTP=  S«RT(KFACT*SQRT(UVMAG) ) 

RATI02=RATI024BTC/BTP 

RAT2=BTC/BTP 

IF(IFLAG.EQ.I)  V AR 2 =V AR 2* < R AT2-C 2 ) ** 2 
C  WRITE (6 »9G1 1 ) 110  ,  J 1 0 , R  AT  2  ,  B  ZC , B  Z  P  *  I T (3»I10»J10) » I T (4,I10,J1Q) « 

C  5  IT(5,I10»J10),IT(6,I10,J1C) 

JCOUNT=JCOUNT ♦ 1 
1  Cl  2  CONTINUE 
1017  CONTINUE 
1  Cl S  CONTINUE 

C  1  P  =  RATI«1/<IC0UNT*1 .) 

C2P=RATI02/(JC0UNT*1.) 

WR  I  TE (6  *1 01 3 )  KFACT,  Cl  ,C 1 P , C 2 , C 2P , I K N  OWN 
1C17  FORMAT (3X,  6H  K  =  ,F5.2,13H  C1(GIVEN)=  ,F6.0,13H  C1(CALC)= 

JFO.Ofl^H  C2(GIVEN)=  ,F6.0,13H  C2(CALC)=  ,F6.C,10H  IKN0WN=  ,15) 
I  F  (IKNOWN . EG . 1 )  GO  TO  1014 

ri=cip 

C<=C2P 

1  Cl  4  CONTINUE 

IF(IFLAG.EO.I)  VARl=SQRT (  V  AR  1 ) /ICOUNT 
If  (IFLAG.E0.1)  VAR2  =  SQRT(VAR2)/JCOUNT 


J 


1  CM  6 

C 

C 

1020 

c 

c 


c 

c 

1030 

C 

c 


c 

c 

1040 

c 


c 

c 


c 


I FLAG=I FLAG+1 
IF (IFLAG.LT. 2)  GO  TO  1015 

WRITE(6,1016)  VAR1.VAR2 

FORMAT (2E1C.3) 

GO  TO  1019 

*  *  *  *  BLOCK  1020 


CONTINUE 


(BLK  140) 


CALCULATE  DIRECTION,  DIRECTION  V  AR  I*AN C E 
PHI=ATAN(BY/BX) 

CHI*-TAN(2.*PHI> 

PEN0M=(I34/I56)*(I56P/I34P) 

DFACT=.5/(1.  ♦  DENOM*DENOM) 

DPD3*  2.*DFACT  *  I T ( 4 , I , J ) *  I  5 6P/ ( 1 3 4P *  1 34P *  1 5 6  ) 
DPD4*  2.*DFACT  *  I T ( 3 ,  I  , J > *  I 56P / ( 1 3 4P* I  34 P *  1  5 6  ) 
DPD5=  2 »*DFACT  *  I T <6 , I , J ) *  1 34  /(l56  *1  56  *I34P> 
DPD6=  2.*DF ACT  *  I T ( 5 , I , J ) *  1 34  /(I56  *156  *I34P) 
SIGPHI  =  SIG32*DPD3*DPD3 
S  +  SIG42*DPD4*DPD4 

$  ♦  SI652*CPD5*DPD5 

1  +  SIG62*DPD6*DPD6 

SIGPHI=SGRT (SIGPHI ) 

IF (SIGPHI .GT.SPHIMX)  S I GPH I =SPHI MX 
GO  TO  1029 

****  BLOCK  1030  **** 


CONTINUE 


(BLK  140 


CALCULATE  B  SIGMAS 

TA*  C2*C2*KFACT*(U*U*SIGU*V*V*SIGV>/(4. *UVMAG*SGRT (UVMAG)  ) 

TB=  (.5/(1.  +  CHI*CHI))**2.*((  S IG3  2*  S IG42  )  /  (  1 34  *  134  ) 

S  ♦(SIG52-»SIG62>/(I56*I56>  ) 

TB*1  . 

TC=SIGPHI*SIGPH1*BT*BT 

SIGBX2  =  TA  *  C0S(PHI)**2  ♦  TC *TB * S I N  (PH  I )  *  *  2 
SIGBY2  -  TA*SIN(PHI)**2  ♦  T C *TB* COS ( PH  I > ** 2 
SIGBZ2* 

$4.*C1*C1*KFACT**2  *  ( S I G 1 2  *  I T (2 , I , J ) **2  ♦  S I G22* I T  ( 1  , 1  ,  J ) ** 2 > / 

SI12P**2/I12P**2 
GO  TO  1039 

**♦*  BLOCK  1040  *  *  *  * 


CONTINUE 

CALCULATE  PARAMETER  SIGM\S,  CHECK  FOR  LARGE  TERMS 
CHECK*. 5 

ST*  BZ*BZ*(DJ*DJ*SIGBX2  ♦  0 1  * D I  * S IGB Y 2  ) 

1  ♦  SIG9Z2*  ( D  J  *  B  X  -  DI*BY)**2 
ST*  D  J  *  D  J  *  (PZ  *9Z  *S I G9X  2  ♦  BX *B X  * S  I GB Z 2 ) 

J  4  D I *D I  *  (B7*BZ*SIGBY2  ♦  B Y *8 Y * S I GBZ 2  ) 

SE*  (BX*DI*BY*DJ)**2*SIGBZ2  ♦  B Z *B Z *  (  S IGB X 2*0  I  * D I *S IGP Y 2  * D J * D J ) 
SFX*  PZ*BZ*SIGBX2  ♦  BX*BX*SIGBZ2 
SFY*  BZ*67*SIGBY2  ♦  BY*BY*SIGBZ2 

SFZ*  4. * (6Z*0Z*SIGBZ2  ♦  BX*BX*SIGBX2  ♦  B Y « P Y * S  I  GB Y 2 ) 


C 


IFd.LT.ICHK.AND.J.LT.JCHK)  GO  TO  1,041 


CKKT  =  SIGT(lO.JO)  *  C  H  E  C  1C 
CHKE  =  SIGE  (  1  0  ,  J  0  )  ‘CHECK 
CHKX  =  SIGFX  ‘CHECK 

CHKY*  SIGFY  ‘CHECK 

CHKZ  =  SIGF7  ‘CHECK 

IFCCHKT.GT.ST  .AND.  CHKE.GT.SE  .AND.  CHK  X.GT.SFX  .AND. 

$  CHKY.5T.SFY  .AND.  CHKZ.GT.SFZ  >  GO  TO  1041 
WRITE (6,1042)1, J,  ST,SE,SFX,SFY,SFZ,BX,BY,BZ 

FORMAT  (215  »36H  FOLLOWING  ARE  T , E  ,  F X  ,  FY  ,  FZ , VAR  I ANC E S , 5E 1 0 . 2 , 3 F8 . 0) 
WRITE (6 ,1043)  SIGT ( I0,J0) , S I G  E ( 1 0 , J  0 ) , S I G  F  X , SIGFY ,SIGFZ 
FORMAT (46H  AND  RUNNING  TOTALS  ARE  ,5E1C.2) 

CONTINUE 

SIGT  (10, J0)=SIGT(I0, JO)  ♦  ST 
SIGE(IC,J0)=SIGE(I0,J0)  ♦  SE 
SIGFX=  SIGFX  ♦  SFX 
S IGF Y=  SIGFY  ♦  SFY 
SIGFZ*  SIGFZ  ♦  SFZ 
GO  TO  1049 
BLOCK  1050  $$$$ 

CONTINUE 

PHI1=ATAN(>U/V) 

IF(V.LT.O.)  PH 1 1 =PH 1 1  +P I 
PH  1 1  =PH 11/ (2.*DT0R) 

°H I =PH I / DTOR 

WRITE(6,1051>  PHI,PHI1,I,J 
PHI  -PH  I  ‘ DTOR 
F0RNAT(2F1C.3,2I5) 

60  TO  1 C59 
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FIGURE  CAPTIONS 

Figure  1.  Distribution  of  polarization  measurements  for  weak  field 
picture  elements  of  solar  active  region  2665,  used  for 
determining  observational  biases  and  uncertainties.  Normalized 
frequency  of  occurrence  is  plotted  against  (a-b)  two  components 
of  linear  polarization,  (c)  circular  and  (d)  total  linear 
polarization.  If  no  biases  have  been  introduced,  distributions 
should  center  and  peak  at  zero.  Small  biases  in  transverse 
measurements  probably  are  induced  by  griding  optics. 

Figure  2.  Distribution  of  polarization  measurements (See  Figure  1)  after 
minor  correction  for  biases.  Spread  in  distribution  (Op  and 
Opg)  gives  random  uncertainties  used  in  estimating  noiseVef fects. 

Figure  3a-b.  Circular  Crosstalk  Determination.  The  measured  circular  polar¬ 
ization  is  compared  with  the  two  components  of  the  linear  polar¬ 
ization,  with  each  plotted  point  representing  a  picture  element. 
Only  the  weak  fields  are  shown  and,  therefore,  the  relationship 
of  linear  to  circular  polarization  should  be  random.  The  weak 
but  systematic  correlation  (seen  here)  between  the  magnitude 
of  the  circular  and  linear  polarization  measurements  implies 
circular  crosstalk  into  the  linear  polarization  measurements. 

In  the  case  shown,  the  deviation  of  best  fit  straight  lines  from 
horizontal  approximates  the  degree  (percentage)  of  crosstalk 
observed.  Although  an  infrequent  problem,  crosstalk  testing  and 
correction  procedures  are  routinely  incorporated  into  the  MSFC 
data  analysis  programs. 

Figure  4a-b.  Circular  Crosstalk  Correction.  Same  as  Figure  3,  after  the 

estimated  correction  for  percentage  of  circular  crosstalk  has 
been  applied  to  the  data.  Linear  polarization  data  now  appear 
free  of  circular  crosstalk  influence.  (Note  that  the  best  fit 
for  the  zero  values  approximates  a  horizontal  line.) 

Figure  5.  Contours  of  circular  polariztion  measurements  for  AR2372  during 
the  time  of  spot  motions  and  greatest  flare  activity.  In  all 
figures,  solid  lines  depict  positive  magnetic  polarity, 
dashed  lines  are  negative,  and  orientations  are  the  same. 

Contours  levels  are  +  0.01,  1.0,  2.5,  and  multiples  of  +  5% 
polarization. 

Figure  6.  Measured  direction  of  horizontal  force  (0)  for  AR2372.  Greatest 
spot  migration  occurred  after  data  points  on  5  April  and  before 
point  on  7  April.  Direction  of  proper  motion  of  the  migrating 
spot  (from  white  light  and  magnetograms,  Figure  5)  agree  well 
with  data  points  on  6  April. 
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Figure  7  Vector  Magnetographs  of  AR2665.  The  line-of-sight  component 
(B^)  is  represented  by  contours  of  field  strength,  solid  for 
positive  values  and  dashed  for  negative  values.  The  transverse 
component  (B^,)  is  represented  by  azimuthal  vectors  where  length 
is  proportional  to  field  strength.  The  minimum  length  repre¬ 
sents  300  Gauss  and  the  maximum  length  represents  1000  Gauss.  Note 
the  changes  occurring  in  the  vicinity  of  the  line  of  B^  polarity 
reversal.  Throughout  the  period  under  study,  there  is  weakening 

of  both  and  the  components  of  the  field.  Continuing  sim¬ 

plification  is  also  shown  by  the  increasing  orthogonal  nature 
(less  shear)  of  the  transverse  component. 

Figure  8.  Histogram  of  X-ray  Flares  for  AR2665.  Logrithmic  plot  with 
Cl  =  10"3  ergs  cm-2  s'"1.  Ml  =  10  x  Cl.  Note  the  decline  in 

activity,  with  only  minor  flares  after  the  12th. 

Figure  9.  Region  History.  The  number  of  spots  is  indicated  by  the  dark 
triangle  in  (a)  and  by  the  asterisks  in  (b).  The  solid  lines 
represent  sunspot  area  evolution.  The  flare  histograms  at  the 
bottom  of  the  plots  indicate  the  number  and  magnitude  of  the 
x-ray  flares;  open  histograms  for  class  C,  lined  for  class  M. 

In  its  first  disk  transit  (a),  the  region  was  very  small  and 
non-productive  of  flares.  As  AR2665  (b),  early  activity  had 
declined  by  mid  transit  and  the  region  exhibited  a  gradual 
decline  in  sunspots  and  sunspot  area. 

Figure  10.  Physical  parameters  (Equations  2-12,  2-11  and  2-13)  for  AR2665, 
determined  using  Lorentz  equations  and  vector  mganetograms 
obtained  as  a  function  of  time  and  observing  wavelength.  Solid 
lines  connect  120  mA  observations,  while  simple  60  and  90  mX 
observations  are  also  shown.  Error  bars  are  one  sigma  values 
determined  as  the  square  root  of  the  sum  over  pixels  of  Equations 
4-1,  4-3,  and  4-2  respectively. 


PTS.  VS  PV  (1A1B)  PTS.  VS  PQ  (SORT  (2AB**2+3AB**2)-PQ0) 


PTS.  VS  PV  (1A1B)  PTS.  VS  PQ  (SQRT  { 2AB* * 2+3AB*  *2) -POO) 


15/1359  U.T.  SEPT 


Figure 


15/1359  U.T.  SEPT 


Figure 


36 


Figure  8 
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TABLE  I. 
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MSFC  Magnetograph  Optical  Configuration  For  Measuring  Circular 
And  Linear  Polarizations 
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